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1. Introduction RES can improve the variety of supply in energy markets, 

reducing thus the dependence on oil; they can also help to 

The global interest in dealing with the energy crisis and safeguard sustainable energy sources on a long-term basis, and 

securing the energy supply has led to the promotion of renewable contribute to a reduction of the local and global environmental 

energy sources (RES), in various sectors of the economy [1]. impact on regional sustainable development. RES account for 17% 


of the global primary energy production, mainly through large 
hydroelectric installations and the use of traditional types of forest 
—— biomass and agricultural residues in developing countries [2]. 
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disquieting increase of carbon dioxide has been observed in the 
environment, while at the same time, its deposits are constantly 
decreasing due to the extensive destruction and degradation of 
natural ecosystems. This fact also contributes to the atmospheric 
pollution of the environment and intensifies climate change [3]. 

The use of biomass as a low cost form of energy could be 
further expanded. However, the absence of a sound energy policy 
aggravates the market situation and puts forward significant 
economic obstacles, which render biomass non-competitive vis- 
a-vis other renewable energy forms and fossil fuels. There would 
be no economic difficulties for the promotion of biomass, if a 
strong policy was put into place that would promote its use, and 
differentiate it from the remaining RES. What is required is the 
implementation of a rational policy that would promote biomass 
use on the global market [4]. 


2. RES and the use of fuelwood 


Fuelwood is one of the main types of biomass and its contribu- 
tion to the energy requirements of a country is significant [5,6]. 

The energy policies of various countries differ depending on 
the needs of each country and their level of dependence on fossil 
fuel imports. Low-income families cannot afford the costs that the 
use of fuel and the purchase of new technological products entail. 
Therefore, an alternative energy policy based on the reduced use 
of fuelwood cannot be promoted in such countries [7]. 

The energy policy objectives been defined by developed coun- 
tries favour the use of RES only for the production of electricity 
and biofuel. Minimum efforts have been made to promote the use 
of fuelwood. The only exception is the widespread use of fuelwood 
in rural areas. Only in recent years, have efforts been observed by 
state authorities to provide the necessary funding for its promo- 
tion. However, given the fact that the relevant funding is minimal, 
no major tax incentives have been instituted [8]. 

The existence of a strategic plan to support renewable energy 
sources would render fuelwood more competitive in the market- 
place, and increase its level of usage compared to other heating 
sources; nevertheless, this would only be feasible if the price of 
fuelwood were lower compared to fossil fuel [8]. 

Fuelwood is a product that is mainly known for its use as a 
heating source. In many places, beyond its essential part in domestic 
consumption, fuelwood is used in a broad range of small businesses 
(bakeries, restaurants etc) at a ratio of 20% of the total consumption 
[9]. Fuelwood usage depends on the quality and type of wood, its 
availability in the natural environment and the social characteristics 
of the populations using it [10]. As an energy source, it has the ability 
to provide a significant percentage of the energy required by low 
income rural families and to ensure a low production cost for the 
urban energy as well. Whether this goal is achieved or not, depends 
directly on the support of local authorities [1]. 

The use of energy from biomass and specifically from fuelwood 
is even today very closely related to the people's quality of life in 
the rural areas of Africa and Asia. Even in countries such as South 
Africa, in which an extensive electrification programme has been 
implemented, wood is still the most widespread type of fuel used 
for cooking [11]. 

In countries such as Bangladesh, and more specifically in its 
rural regions, biomass is the primary energy source for house- 
holds. In order to cover their energy requirements, household 
owners mainly collect bamboo, branches, cow dung, fuelwood, 
rice husk, leaves, twigs and straw [12]. 

In addition, wood is still the type of fuel chosen by most urban 
households in the capital of Burkina Faso. The demand for energy 
from wood is directly linked to household income. The percentage 
of fuelwood usage decreases as household income rises [13]. 


In an area of Uganda, it was found that there is almost full 
dependence on the use of fuelwood for household cooking and for 
small industries. Most of the fuelwood for household use is 
collected by women and is mainly dead wood, resulting thus in 
a limited environmental impact on forest ecosystems [14]. 

In Greece the energy consumption is expected to rise significantly 
over the next few years, since it is considered essential for the 
country's economic growth, if it is to cover the gap that separates it 
from the developed economies of the European Union [15]. 

The country however is dependent to a large extent on the 
import of oil and natural gas. The increased energy requirements 
of Greece cause significant problems to the balance of payments, 
with a further impact on the other macro-economic figures of its 
economy [15]. There is therefore a need to further intensify the 
exploitation of domestic energy sources, and particularly RES. 

The share of RES in the national energy balance is around 5.6% 
of the total gross domestic consumption and 17.7% of the domestic 
primary energy production. In 2008, the production of primary 
energy from RES was 1.8 Mtoe, compared to 1.2 Mtoe in the early 
1990s. Of these, 600 ktoe approximately are attributed to the use 
of biomass in households and 264 ktoe to the use of biomass for 
industrial needs [16]. 

Although the overall contribution of RES to the gross domestic 
energy consumption is relatively low (5.5-6.5%) the contribution 
of RES other than household consumed biomass and large hydro- 
electric plants is presenting a steady increase due to measures 
targeted to its economic support [16]. 

In Greece, the primary energy form of biomass is the wood 
produced by Greek forests, which is mainly used as fuel by 
households. Agricultural crop residues constitute also an impor- 
tant and promising biomass source (branches, straw etc). 

Greek forests and other forest areas produce a range of 
products and services [17,18]. Fuelwood is one of the primary 
products of Greek forests, since it accounts for over 65% of the total 
wood production in recent years [19,20]. The production of fuel- 
wood mainly comes from oak forests, which cover a large part of 
Greece [21]. In recent years, there has been a shift in European 
Union policies regarding the development of rural areas (agro- 
environmental measures, LEADER+ etc), which has opened up 
new opportunities for the exploitation of agricultural land by its 
owners [22,23]. Thus, the plantations of forest species that have 
been created in abandoned or low productive agricultural lands 
within the framework of EU policies seem to be in position to 
contribute to the increase in the produced quantities of fuelwood 
over the next few years [24]. 

Total wood production has significantly decreased in recent 
decades (from 1.7 million m? in 1990 to 0.96 million m? in 2007 
[25]). This reduction had largely affected fuelwood compared to 
industrial wood [26], due to the much better prices paid to loggers 
for the sale of industrial wood. However, the economic crisis of 
recent years lowered the demand of industrial wood and increased 
the demand for fuelwood; nowadays significant quantities of even 
industrial wood are being sold to households as fuelwood. 

In the past, the consumption of fuelwood was very high in rural 
areas, particularly those located in mountainous and semi- 
mountainous parts of the country. Its consumption had gradually 
decreased following the reduction of rural population (due to 
urbanism), the improved standards of living and the substitution 
of fuelwood with other fuel, primarily oil. However, the oil crises 
and the use of fireplaces in urban residences and country houses, 
had contributed somewhat to an increased demand for fuelwood 
[26,27]. Today, due to the economic crisis affecting the country, the 
consumption of fuelwood has increased in a number of house- 
holds, mainly those situated in rural areas [28]. 

[Forests in Greece, according to the latest census, cover 
33,58,000ha or 25.4% of the country's total area. Approximately 
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2/3 or 65.4% belong to the state and the remaining 34.6% are not 
state-owned, and belong to private entities, Local Authorities, 
monasteries and other public benefit organizations [29]. The body 
responsible for the exploitation of forests is their owner as a rule. 
Therefore, the body responsible for the exploitation of public 
forests is the Greek state. Until the end of 1986, two systems for 
the exploitation of public forests were used in Greek forestry. One 
was the self-supervision system, within the framework of State 
Forest Exploitation (KED, from the Greek acronyms) by the local 
Forest Services [18], according to article 137 of L.D. 86/69, and the 
other was through the leasing of the forest's production (timber 
cut) by forest cooperatives, according to article 134 of L.D. 86/69 or 
following an auction among those interested in leasing, according 
to article 120 of L.D. 86/69 [18,30]. 

Through the first system, the Public (Forest Service) carries out 
the exploitation of forests by assigning their logging to the Forest 
Cooperatives, who are paid according to the timber produced. 
According to the second system, the Public leases the production 
to Forest Cooperatives, who pay a lease depending on the volume 
of the timber expected to be produced, and then sell the wood 
products themselves. 

With article 74 of L. 1541/85 and of P.D 126/86, a new system 
of exploitation (the third) was introduced for the exploitation of 
public forests, according to which the exploitation of forests is 
assigned to Agricultural Forest Cooperatives, who pay a percentage 
of the selling price of forest products to the Central Fund for 
Agriculture, Livestock and Forests and to the Local Authorities, 
within whose administrative boundaries the specific forest, forest 
part or stand is situated [18,30]. 

In this paper a production—allocation model is presented for the 
optimal management of fuelwood in Greece. The need of a 
production management model that would incorporate both 
environmental and economic criteria is more imperative now, 
due the rapid increase in fuelwood demand. In a fuelwood supply 
chain, the “factory” is considered to be the forest. Through the 
channels of the supply chain (AFCs and fuelwood merchants) the 
product arrives at the customer (households, industries etc). 
Fuelwood production and demand are linearly dependent, so 
when there is an increase in demand, the production should be 
also increased, so as to cover the demand. Increased fuelwood 
production affects the sustainability of the natural resource, 
leading to irreversible situations. The proposed model, which is 
an extension of [31], incorporates the restrictions posed by 
legislation (P.D. 126/86) on the maximum produced fuelwood 
quantities, minimizing the overall cost (production, transporta- 
tion, and fixed installation cost of the warehouses). 

The rest of the paper is organized as follows: In Section 3 the 
problem statement and the proposed mathematical formulation 
are presented and discussed and in Section 4 expected results and 
further research on the proposed model are presented. 


3. Mathematical formulation 
3.1. The fuelwood supply chain 


The fuelwood supply chain (FSC) is quite similar to this of a 
common product but with a few changes. Examining the supply 
chain in a general form several nodes are considered. A FSC is very 
important to rural development and can be characterized as a 
form of local bioenergy system (LBES) [32]. At the beginning of the 
supply chain there are the production plants, the units that 
combine all the raw material for the manufacture of the product. 
The next step after production is to transport the product to an 
intermediate node that is used for storing and for distributing the 
product to the customers, depending on the demand. In the final 


stage of the supply chain there are the customers whose demand 
must be satisfied [33]. In order for the companies to achieve the 
best level of customers’ service they try to improve not only the 
design of the supply chain (SCD) but also the supply chain 
operations (SCO) conducted amongst the different nodes. In the 
present case the above formulation can be used with a few 
changes. Each company depending on the nature of the product 
and on the structure of its supply chain considers one of the 
following supply chain models: the pull or the push systems [34]. 
The first systems are also known as centralized systems, as central 
decisions are made about the produced quantity. In a centralized 
system the production is planned without taking into considera- 
tion the demanded quantity of the product. This production policy 
may lead to overbalance of the demanded quantity. The latter 
production systems, also called decentralized, adjust the produced 
quantity according to the information brought from the customers' 
demand. These systems adopt built to order (BTO) policies and can 
be well fitted to computers or car production systems. 

The first step in setting up a FSC is to replace the production 
plants with forests, which are the production “facilities” for the 
chain. Due to the fact that a forest has a limited produced quantity 
and a limited time horizon for the cutting, a pull system for the 
FSC can be adopted. Each forest is divided, for managerial reasons, 
into forest compartments. The compartments to be exploited each 
year, according to the management plan, are assigned to Agricul- 
tural Forest Cooperatives (from now on AFCs) following the 
regulations of laws mentioned previously. AFCs are located into 
villages inside or in the vicinity of the forest. The wood of felled 
trees is not formulated into fuelwood at the forest compartments, 
but it remains in the form of relative long logs for cost effective 
transportation. The produced quantity in this form is directly sold 
by the AFCs to fuelwood merchants, which in turn cut it into small 
logs appropriate for consumption and sell it to customers. Facilities 
of the AFCs and merchants are considered already installed with a 
given capacity. The amount of fuelwood demanded from the AFCs 
and the merchants is complementary. Quantities of fuelwood are 
considered to be delivered to customers' zones which represent 
large areas, as villages or cities. The demand for fuelwood is 
estimated in advance, based on statistical data from previous 
years. The examined FSC is presented in Fig. 1. 

Main objective of the study is the minimization of the total cost 
of the processes occurring in the FSC. As the proposed model is 
demand driven (when demand increases the supply should rise at 
appropriate levels to cover the demand) an increase would lead to 
overall increase of the quantities produced and therefore to the 
overall cost. Thus, in order to avoid forest degradation and over 
logging that rapid changes in fuelwood demand may create, the 
fuelwood quantities produced are assumed to be upper and lower 
bounded. Due to the large distance between merchants and 
customers and the specific conditions for fuelwood maintenance 
(to guarantee low moisture and easy and effective combustion) 
intermediate warehousing facilities located at merchants’ sites 
are also assumed. For this purpose a mixed integer linear 


et Merchant }———> ‘Customer 


m 


Forest 
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Fig. 1. The examined FSC. 
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programming (MILP) model is examined providing information 
about the following processes throughout the planning horizon of 
the considered FSC: 


© The fuelwood volume produced. 

è The levels of wood transported from AFCs to fuelwood 
merchants. 

The levels of wood transported from merchants to customers. 
The size of warehouse, where fuelwood is stored. 

Location and number of warehousing facilities. 

The total cost of all the processes of the examined FCS. 


With the proposed MILP model, decisions about the produc- 
tion, distribution, storage of just fuelwood quantities are proposed, 
with the use of Lagrangean relaxation algorithm. The model is 
oriented towards the Greek P.D 126/86 for forest production 
management (such model has not yet been proposed for produc- 
tion management in Greece). Moreover, due to the fact that any 
changes in fuelwood demand are incorporated in the proposed 
model, an overall picture about the cost of the operations 
conducted and the levels of quantities produced and distributed 
throughout the supply chain are also derived. Through the use of 
algebraic formulation the previous structure can have multiple 
echelons at each stage of the examined FSC. This provides to the 
proposed model a more realistic approach to the real world 
problems (Tables 1-4). 


Table 1 
Indices and sets. 


peP Forest compartment 
feF Forest 

seS AFC 

meM Merchant 

ceCU Customer 

ieSC Scenario 


Table 2 
Continuous and decision (binary) variables. 


Qpfs Quantity of fuelwood logged at compartment p of forest f by AFCs (m?) 

P;; Cumulative quantity of fuelwood logged of all compartments of forest f by 
AFC s (m?) 

Fm Quantity of fuelwood processed at merchant's node m (m°) 

Qsm Quantity of fuelwood sent from AFC s to merchant m (m°) 

Qme Quantity of fuelwood sent from merchant m to customer c (m?) 

Wm Capacity of warehouse installed at merchant's location m (m°) 

Ym 1 if warehouse m will be installed in merchant's m position, 0 otherwise. 

Xmc 1 if the connection between merchant m and customer c exists, 


3.2. Forest production constraints 


At the initial stage, there is the forest that is the productive unit 
of the supply chain. The land use should be subjected to strict 
constraints regarding sustainable management and conservation 
of the environment [35]. Fuelwood productive forests are regu- 
lated even-aged forests, which are managed under a sustainable 
framework to ensure their long term productive, protective and 
environmental services. Therefore, restrictions to the harvesting 
quantities from each forest stand are set by the management plan 
in order to both maintain the productive ability of the forest and 
sustain its non-timber functions. All quantities mentioned below 
are considered to be converted from stacked cubic meters (st m?) 
to true cubic meters (m?) following the next conversion: 1st 
m?=0.6 m’. 

Let Qprs be the continuous variable that represents the amount 
of wood harvested from compartment p of the forest f from AFC s. 
The quantity harvested from each partition is bounded between an 


upper ( pfs) and a lower (Q£,.) amount. This leads to the following 


: pfs 
constraints: 
Qoys<Qprs> YP. f, S (1) 
Qoys>Qhys> VD.F.S (2) 


Consequently the quantity harvested from all the forest's 
compartments would yield the production of fuelwood of forest 
f. This amount is represented by the continuous variable Pps, and 
models the amount that will be transported through the channels 
of the FSC and will be finally delivered to the last node, the 
customer. Hence the following constraint is introduced: 


> Qpfs = Pfs, vf,s (3) 


Table 4 
Cost parameters. 


cout Cost of cutting trees. (€/m?) 

Mov Cost of moving from forest's compartment to the forest road (€/m?) 
MUN Fixed tax placed at AFCs according to P.D. 126/86 (€/m°) 

GOV 

PURCH Purchase value of fuelwood quantity by merchants (€/m*) 


CR ‘i Transportation cost of fuelwood quantity from AFC to merchant 


(€/m?) 
cre Processing cost of fuelwood at merchant's node (€/m°) 
cr. é Transportation cost of fuelwood from merchant to customer (€/m?) 
cn Fixed installation cost of a warehouse to merchant's location (€) 


TCFSC Total cost of the examined fuelwood supply chain (FSC) (€) 


0 otherwise. ETCFSC Expected total cost of the examined fuelwood supply chain (FSC) (€) 
Table 3 
Capacity and demand parameters. 
Qu. Upper bound of quantity of wood logged at compartment p of forest f by AFC s according to P.D 126/86 (m°) 
Qt a Lower bound of quantity of wood logged at compartment p of forest f by AFC s 
P% Upper bound of cumulative quantity of fuelwood logged of all compartments of forest f by AFC s (m°) 
Pr, Lower bound of cumulative quantity of fuelwood logged of all compartments of forest f by AFC s (m?) 
Qu. Upper bound of quantity of fuelwood transported from AFC s to merchant m (m°) 
Qhn Lower bound of quantity of fuelwood transported from AFC s to merchant m (m?) 
Qu. Upper bound of quantity of fuelwood transported of merchant m to customer c (m?) 
Qh, E Lower bound of quantity of fuelwood transported of merchant m to customer c (m°) 
De Demand for quantity of fuelwood at customer c (m°) 
wY, Upper bound for the capacity of warehouse installed at merchant's location m (m°) 
wt Lower bound for the capacity of warehouse installed at merchant's location m (m?) 


Mi Occurring probability of scenario i. 
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Produced quantity is also bounded between an upper (PH) and 
a lower (Pi) amount. The previous bounds can be derived from 
the summation of the upper and lower bounds of the compart- 
ments of each forest as shown in Eqs. (4) and (5): 


Pr = È Qos Y.S (4) 
P 
Pi = x Qis Ys (5) 


3.3. Fuelwood balance constraints 


The amounts of fuelwood produced at the forest are transferred 
though the channels of the FSC to the end customer. At the 
beginning of the supply chain the produced volume of fuelwood 
from forest f by AFC s, should be equal to the quantity of 
fuelwood's volume transported by AFC s to merchant m. Let Qsm 
be the quantity of fuelwood transferred from AFC s to merchant 
[36]. The previous statement is modeled by the next equation: 


Pts = È Osm: vs (6) 


From the node of AFC the flows of fuelwood are transported to 
the customers' zones. So the quantity delivered from AFC s to 
merchant m should be equal to the quantity delivered from 
merchant m to customers' zone c. Let Qmc be the continuous 
variable that represents the quantity of fuelwood delivered from 
AFC s to customer zone c then the next constraint is introduced: 


2 Qsm = x Qmo» vm (7) 


As the transportation of quantities of fuelwood from AFCs to 
merchants is being done with trucks of finite capacity, the 
delivered quantities cannot exceed the maximum capacity of the 
trucks ( Ti: Moreover the delivery of quantity of fuelwood can be 
done if the amount exceeds a minimum certain level (Q{,,). The 
previous lead to the next constraints: 


Qsm<Qy,, VS, m (8) 


O20 as YS, m (9) 


Fuelwood quantities that reach the merchant's node are appro- 
priately processed in order to ease the customers. The quantity of 
fuelwood processed at merchant's node (Fm) can be computed 
with the following constraint: 


Fm=YQsm, vm (10) 


The quantity of fuelwood transported from merchants to 
customers’ zones should be equal to the demanded quantity of 
fuelwood at customer zone c (D+). The previous can be modeled 
through the following constraint. 


È Qc = De, vc (11) 


In merchants' node, the number of warehousing facilities is 
determined through MILP solution. Let Ym be a decision (binary) 
variable, that take a value of 1 if warehousing facility (Wm) will be 
installed in merchant's m location and 0 otherwise. The possible 
installed facilities have finite capacity and have an upper and a 
lower bound. The previous proposition is modeled with the 
following constraints: 


WinsWY x Ym, Ym (12) 


Wm>W, x Ym, Ym (13) 


As in constraints (8) and (9), the quantities delivered from 
merchants to customers should be bounded between a lower 


(QE) and an upper (QUe) bound concerning the minimum amount 
that can be delivered and the maximum capacity of the trucks that 
are used for the transportation, respectively, and are delivered 
only if the connection between warehouse m and the correspond- 
ing customer, exist. If Xmc is a decision variable that represents the 
connection between merchant's warehouse m and customer c, the 
next constraints are introduced: 


Qme<SQme x Xme, YM, c (14) 
Qinc= Qn x Xme, YM, C (15) 
Xmc<Ym, VM, C (16) 
VXme = 1, Ve (17) 
m 


Inequality (16) implies that a connection between the ware- 
house placed located at a merchant's site would exist if-f ware- 
house m is installed. Equality (17) implies that each customer c is 
served by exactly one warehouse m (single sourcing constraint). 


3.4. Non-negativity constraints 


As all the continuous variables introduced in the previous 
paragraphs represent quantities of fuelwood, by solving the linear 
programming problem the resulting values of these continuous 
variables cannot be negative. Therefore the next non-negativity 
constraints are introduced: 


Qyys20, Vp.f.s (18) 
Prs>0, Vf,s (19) 
Qsm=0, Ys, m (20) 
Qmc20, Ym, c (21) 
Fm>0, Ym 

Wm20, Ym (22) 


3.5. Objective function 


As main objective of the proposed mathematical model is to 
minimize the total cost of the operations of FSC, in this section the 
terms of the objective function will be analyzed. 

The forest is considered the first link of the examined FSC and 
different compartments of the forest are exploited by AFCs. Let 
CT represent the cost of cutting the trees and CM9” of moving the 
trees from forest's compartment to the forest road. The overall cost 
of the previous operations is modeled with the next term: 


(cour + coy, J IDX Qpfs (23) 


As mentioned above the total amount of fuelwood logged from 
all compartments of a forest is equal to the amount of fuelwood 
that is produced from the forest. According to P.D 126/86, a tax is 
placed on AFCs for the exploitation of the forest. A part of this tax 
is attributed to the municipality where the forest geographically 
belongs to (MUN), another part is attributed to the government 
(GOV). Based on the above, the next term is introduced in the 
objective function: 


(MUN + GOV) x X > Pjs (24) 
Ff Ss 


The next link after AFCs is the merchants who purchase the 
logged quantity from AFCs. Here also a variation of operations is 
done regarding transportation, and cutting the pieces in smaller 
ones. Let PURCH be the purchase cost of logged amount of 
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merchants from AFCs and Cem be the cost of transporting the 
quantity of fuelwood from AFCs to merchants. So the next term is 
introduced in the objective function to model the transportation 
cost: 


(PURCH + Carc.m) X X È Qsm (25) 
s m 


At merchant's node one of the operations that are done is the 
cut of the logged quantity into smaller pieces so that would be 
easy to use in a household fireplace. If CER denotes the processing 
cost of the above described operation done at merchant's node 
then the following term is introduced in the objective function: 


Cre x E Fm (26) 
m 


The final operation at the merchant's node is that of the 
transportation of fuelwood to the customers. Let Che represent 
the transportation cost of fuelwood quantity from merchants to 
customers. Thus the following term is introduced in the objective 
function: 


Cr se LL Qne (27) 


For the installation of a warehouse facility in a certain location, 
a fixed installation cost is required. Let C/V be the fixed installation 
cost of the warehousing facility in merchant's location m. The total 
installation cost is expressed with the following term in the 
objective function: 


ECN x Vin (28) 
m 


The total cost of the examined FSC(TCFSC) integrates all the 
costs of the operations conducted within the chain. Hence the final 
form of the objective function is derived from the summation of 
terms (23)-(28). 


minTCFSC = (CU? + CM) x 2 2 x Qpfs 


+(MUN + GOV) x © ¥ Pjs 
f S 


+> CN x Yin (29) 
m 


3.6. Implemented policy for the use of fuelwood 


Fuelwood is one of the major fuels for heating of households in 
Greece, because of its low price, which is among the lowest found 
in EU countries [37]. Also the rapid drop of the average income per 
capita and the increase of the heating oil price make fuelwood a 
growing market in Greece. Of course this may have negative 
impacts on the management alternatives available for the forest, 
because as the demand increases, the pressure for higher produc- 
tion quantities should increase as well. Examination, therefore, of 
different demand situations cannot only lead to conclusions on the 
size of the market to be developed but also give an estimate of the 
pressure to be exerted to the forest management. 

Consequently, a scenario based linear programming model is 
adopted in order to integrate the expected changes in fuelwood 
demand [36]. Each scenario, which corresponds to an increase or a 
decrease of the fuelwood's demand, is weighted by an occurring 
probability in the objective function and a superscript s is inserted 
at each continuous variable of the model. Let DË be the demand of 
fuelwood under different scenarios, i=1,2,...,SC. Based on the 
previous, the scenario based demand is the following ICI x ISCI 
tableD® = [D D?’D®---DS°], where le! measures the size of set o. 


Each scenario is considered to be complementary, so the corre- 
sponding probabilities should verify the next equation: 


» m= 1 (30) 
1 
As continuous variables that concern amount of fuelwood are 


subjected to the variations of fuelwood demand, a superscript is 
introduced. So all the continuous variables of the examined model 


are the following:Q° , pO, ® and QË.. Based on the previous 
constraint (3) will become: 
YOU =P) sift (31) 


P 


Of course the conception of the different scenarios does not 
affect only the continuous variables of the model, but also the 
parameters. Regarding productivity upper bounds, following the 
management plan of a forest, the quantity of fuelwood exploited 
by AFCs is given and cannot be exceeded if demand for fuelwood 
increases. Thus the upper bound parameters for any given scenario 
i are considered constant and equal to: 


Qe = Dares WPS Si (32) 


po = PH, Yf,s,i (33) 


Regarding the amounts transported through all nodes of the 
examined FSC also the upper bounds are not subjected to any 
changes. This is happening due to the fact that the quantity of 
fuelwood transported is done by trucks with finite space thus the 
upper bounds cannot change as demand changes. The previous 
will lead to the following equations: 


YO Qu. vs,m,i (34) 


me = Qmo WM, c,i (35) 


As costs of the examined FSC are fixed, the objective function 
measures the expected minimum cost, as each scenario is 
weighted with a corresponding probability. Thus the new objec- 
tive function represents the expected total cost of the examined 
FSC and can be written in the following form: 


minE(TCFSC) = X ro x (C + OM) YY or. 
i Pf Ss 


+(MUN + GOV) x X > PÈ 
fs 


+(PURCH + Chic .m) x EL QS, 
s m 
+C x Z Fn + Cmc x Z Z Oe 


+E C x Ym (36) 
m 


4. Solution approach 


The fuelwood supply chain, production—allocation model pre- 
sented above, is a MILP model. The convex problem consists of 
objective function (29) and linear constraints (1)-(22). A MILP 
problem can be formed in its general form as follows: 


P:min z=f'y+c'x 
s.t 

Ax>b 

Dy>d 

x>0 

ye{0, 1} 


(37) 
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The Lagrangean relaxation problem (LRp) of (37), also called 
Lagrangean dual, has the following forms: 


LRp : min z = f'y + u'(b—Ax) 
u>0 

S.t 
Dy>d (38) 
x>0 
ye{0, 1} 

Using Lagrangean relaxation technique feasible solutions are 
obtained that lead to further minimization of the objective 


function. 
Therefore the following inequality is valid [38]: 


ZRMILPSZLRp ŠZP (39) 


From (39) Zrvmp represents the value of the objective function 
of the relaxed MILP problem. Dualizing single sourcing constraint 
(17) such that:1-2mXmc, the Lagrangean relaxation problem is the 
following [39]: 


maxLRp =min (TcFsco% +> Wo) x (1-2%ne) ) 
s.t 
XE” 


(40) 


In problem formulation (41), x represents the variables of the 
problem and øp is a set that is formed from constraints (1)-(16). 

The steps of Lagrangean relaxation algorithm [40,41], are 
presented below. 


Solve P(Objective function (24), Constraints(1)-(17)) 
Upper bound : L*, Initial value for u* 

05 =2 

for j=0,1,....,.M 

7 =g) 
Uj+ı = Max{0, uj + Uy} 

if |[Uj1-Ujll < € 

Stop 

end if 

If no progress in more than M' iterations 
Gj41 = 65/2 

else 


| Initialization 


(Step size) 


Subgradient iterations 


Gin = 8 
end if 

j=j+1 
end for 


In Lagrangean relaxation algorithm presented previously, the 
following assumptions regarding the initialization of the algorithm 
are made: 


a. The upper bound provided for this algorithm is assumed to be 
equal to the value of the objective function of the relaxed MILP 
(RMILP) problem, P. 

b. The initial value provided for u* is assumed to be the marginal 
value of constraint (17). 


4.1. Description of the numerical example 


In this section a numerical example of the proposed MILP 
model is exhibited. A medium instance regarding the number of 
compartments, forests, AFCs, merchants and customers is exam- 
ined. For this numerical example it is assumed that the area for 
logging covers two forests. For sustainability reasons, fuelwood is 
produced from two compartments which are considered to have 


Table 5 
Capacity and cost parameters of the numerical 
example. 
Parameter Value ( x 1000 m?) 
u 
a, U[5,000, 10,000] 
L 
Qs U[200, 500] 
Pfs EpQ;rs 
Pis LpQhys 
Qu, U[5,000, 30,000] 
at, U[1,000, 1,500] 
Qu. U[5,000, 30,000] 
at, U[500, 1500] 
De U[10,000, 20,000] 
Cost Value 
CUT 4 Mov 100€/1,000 m? 
MUN + GOV 50€/1,000 m? 
PURCH + CRo m 130€/1,000 m? 
cR 100€/1,000 m? 
cN 100,000€ 
27450000 
27400000 
27350000 
S 
© 27300000 
[2] 
ea 
© 27250000 
z 
° ~ 
æ= 27200000 
ks 
27150000 
27100000 
27050000 
TCFSCLR TCFSCRMLP TCFSCmiip 


Fig. 2. Total FSC cost for different methodologies used. 


the right age for logging. Quantities of fuelwood are transported 
by three AFCs over this area to the four merchants of the example 
so as to cover the demand. The last link of the FSC is the customers 
which are not treated as single entities but rather than groups of 
people in villages of small cities. The proposed model is generic 
and can be applied to any number of nodes of FSC. 

The examined model is in a steady state form (time is 
eliminated from the model). For this reason, cost values, upper 
and lower capacity of transported quantities and fuelwood 
demand, are uniformly distributed and are assumed to represent 
annual mean values. All the aforementioned values regarding cost, 
upper and lower bounds, capacity and demand are presented in 
Table 5. Quantities regarding supply and demand express 1000 m?. 


4.2. Results 


Applying Lagrangean relaxation algorithm to P, based on the 
assumptions made above, the minimum TCSFC cost was obtained 
after 5 iterations. Fig. 2, presents the overall FSC cost (TCFSC) with 
the use of Lagrangean relaxation algorithm, solving (37) as an 
MILP and as a relaxed MILP. 

From Fig. 2 inequality (39) is valid as the value of minimum 
cost provided by RMILP, is in this case equal to the value of the 
objective function using Lagrangean relaxation algorithm. Solving 
the primal MILP problem, as it can be seen in Fig. 2, yields a 
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Fig. 3. Solutions of P, solving problem (37) with and without Lagrangean relaxation algorithm. 
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Fig. 4. Solutions of Qsm solving problem (37) with and without Lagrangean relaxation algorithm. 


solution which can be further minimized. From the previously 
methodologies developed, the quantities of fuelwood produced 
and transported throughout the FSC network is of great interest. 

In Fig. 3, the solutions of the problem obtained with the use of 
Lagrangean relaxation regarding the production of fuelwood from 
each compartment by AFCs (Pfs) show a decreasing but smooth 
trend. On the contrary, the solutions derived from the primal 
problem, present large fluctuations which may be more econom- 
ical in some cases but when dealing with natural resources, 
sustainability is not guaranteed. The same stands for solutions of 
variable Qsm that represents the quantities transported from AFC s 
to merchant m, as it can be seen in Fig. 4. Rapid fluctuations are 
presented in MILP solution. 

In Fig. 5, fuelwood quantities delivered from merchant m to 
customer c, are presented. The solutions presented in Lagrangean 


relaxation, show a constant flow of fuelwood from all merchants 
to all customers, unlike the solutions derived from MILP approach. 
The previous is attributed to the fact that in the primal MILP model 
constraint (17) is considered active, while in Lagrangean relaxation 
algorithm this constraint is dualized and added to the objective 
function, as stated in problem (39). 


4.3. Sensitivity analysis 


Due to the nature of the product examined in this work and the 
rapid increase in demand, several distinct cases of demand 
representations (scenarios) are examined. At each demand repre- 
sentation a corresponding probability is assigned. Thus the pro- 
posed cases are weighted into a single objective. In this case, 
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Fig. 5. Solutions of Qmc solving problem (37) with and without Lagrangean relaxation algorithm. 


Table 6 
Scenarios of increase in demand. 


Scenario Demand increase 

1 No increase (De) 

2 Conservative increase in demand by 50% (1.5 x De) 
3 Extreme increase in demand by 100% (2 x De) 


problem (39) is modified as follows: 


E(TCFSC) = X a(i) x TCFSC(x') 


i 


(41) 


In Table 6, the demand representations are presented. The first 
scenario states that fuelwood demand will not increase in the 
future. The second scenario states a quite conservative increase in 
fuelwood demand by 50%. Finally the third scenario states an 
extreme increase in demand by 100%. 

The choice of weight of each demand scenario was made 
towards larger fuelwood demands. Let 2(1), 2(2) and (3) be the 
corresponding probabilities of: (a) the current demand, (b) of the 
conservative increase in fuelwood demand and (c) the extreme 
increase in fuelwood demand, respectively. Three cases of prob- 
ability ranges are proposed for this analysis measuring the impact 
of each demand representation on TCFSC: 


© z(e, | )<(0,0.15], small 
è z(e, -)<(0.15, 0.4], medium 
© z(e, t)<(0.4, 1], large 


The notation z(e, *)used for this analysis states that probability 
corresponding to a demand representation e, is weighted with an 
* impact on objective function. Two distinct set of probabilities are 
proposed. The first one isz(1, 4), 2(2,-) and 2(3, t) which implies 
that extreme increase is weighted more than the conservative 
increase while current demand representation is weighted with a 
small probability. The second one is z(1, 4), z(2, 1) and x(3,-) 
which implies that the conservative increase in demand is 
weighted more than the extreme increase in demand, while 
current situation in fuelwood demand is weighted with only a 


a 
2.90E+07 
285E07| mL 
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Fig. 6. Cases of scenarios weights on objective function using Lagrangean relaxa- 
tion (LR) algorithm (grey line) and solving the MILP model without specific 
algorithm (black line) for: (a) 2(1, +), 2(2,-), 2(3, t), (b) 2, 4), 22, t), 23, -). 


small probability. In the two aforementioned cases, current situa- 
tion is weighted with only a small probability as main aim of this 
analysis is to examine the impact of an increase in demand on 
ETCFSC. Selection of weights is subjected to constraint (30). 

In Fig. 6 a comparative analysis is conducted measuring the 
reduction of the overall cost of the supply chain, solving the 
problem without any algorithm (MILP) and using Lagrangean 
relxation algorithm (LR). As it can be seen in Fig. 6 case (a), 
expected cost has a high value for larger weights placed on 
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extreme demand increase scenario. As weights on extreme 
increase scenario decrease, while conservative scenario increases, 
the expected total cost drops to low levels. In Fig. 6 case (b), there 
is a rapid drop in expected total cost when conservative scenario's 
probability decreases while extreme scenario's probability 
increases. Weighting more on the extreme demand increase 
scenario leads to extreme increases in expected total cost of FSC. 
For both cases of probabilities sets, it can be seen that the 
difference between the two costs derived is constant. The solution 
of the problem with Lagrangean relaxation algorithm yields lower 
ETCFSC for all of the instances of demand representations and 
each scenario. 


5. Conclusions 


Today, the policies of a large number of countries are focusing 
on energy security and efficiency, particularly those whose depen- 
dence on imported energy resources is rapidly increasing. 

The European Union encourages and reinforces the valorization 
of RES, and specifically biomass, for energy production, particu- 
larly when it is compatible with its policies in other sectors, such 
as the common agricultural policy, its environmental policy and 
regional policy. Even today, several countries continue to use wood 
as fuel extensively, for heating and cooking. 

During the last few years, as a result of the economic crisis in 
Greece, a consumer preference shift has been observed towards 
fuelwood. The exploitation of fuelwood contributes to a reduction 
of the country's dependence on imported fuel, improves its trade 
balance and saves foreign currency. Therefore, the sustainable 
management of fuelwood production, the production of high- 
quality products, the minimization of transportation costs and, in 
general, of production costs are some of the basic aims of 
forest farms. 

More specifically, an effort is made in this paper to formulate 
the FSC problem via a MILP model. Lagrangean relaxation algo- 
rithm is applied for the solution of MILP problem leading to great 
minimization of overall cost. As discussed above the model can 
provide information about the fuelwood quantities that are 
necessary so as to satisfy an extremely varying demand. Compar- 
ing the solutions approaches, it can be concluded that when 
applying Lagrangean relaxation algorithm to the presented MILP 
model, the derived quantities present small variance than the 
solutions derived from the solution of the MILP model. Due to the 
fact that demand of fuelwood is not constant, but uncertain, 
scenarios are used to model the increase or the decrease in 
fuelwood demand. Thus if a rapid increase in demand may occur 
this will happen with a corresponding probability. If this prob- 
ability is high then this scenario is more likely to happen, whereas 
the other scenarios are less likely because their probabilities are 
complementary. In this case the model will weight all the 
scenarios according to the occurring probability and the fuelwood 
quantities that will satisfy all demand scenarios will be produced. 
As the primary objective of the model is the minimization of the 
total operational cost, each scenario corresponds to a cost which in 
turn is integrated through probabilities to the overall expected 
cost. For the two sets of probabilities tested, it can be concluded 
that extreme fuelwood demand can eventually increase more the 
expected cost than the conservative increase in demand. 
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